On the critical level-curvature distribution 
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The parametric motion of energy levels for non-interacting electrons at the Anderson localization 
critical point is studied by computing the energy level-curvatures for a quasiperiodic ring with twisted 
boundary conditions. We find a critical distribution which has the universal random matrix theory 
form P(K) ~ l^^l""^ for large level-curvatures \K\ corresponding to quantum diffusion, although 
overall it is close to approximate log-normal statistics corresponding to localization. The obtained 
hybrid distribution resembles the critical distribution of the disordered Anderson model and makes 
a connection to recent experimental data. 
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The Landauer-Buttiker(LB) scattering approachflj 
nowadays is almost exclusively used for computing the 
conductance g for many materials in nanoscience. This 
very appealing formalism although not extremely rigor- 
ous is particularly suitable for small systems, something 
which became obvious following the discovery of conduc- 
tance steps in ballistic point contacts more than a decade 
ago In such a scattering approach one gets the non- 
equilibrium transport for any suffiently realistic system 
if it is connected to two ideal leads. The conductance is 
simply the transmission probability through the sample 
scatterer which can be evaluated at any energy monitored 
by gates, e.g. via Green's function techniques. Although 
g is strictly speaking a property of open systems it can 
be also probed for a closed system from the eigensolu- 
tions of the stationary Schroedinger's equation. This is 
done via Thouless's intuitive definition of conductance 
which involves a measure of the sensitivity of eigenvalues 
to twisted boundary conditions (BC) In the later ap- 
proach, the conductance is related to the so-called level- 
curvatures and is usually expressed by the ratio of the ge- 
ometric mean of the absolute curvature at zero twist over 
the local mean level-spacing. The scattering approach 
for open systems and the eigenvalue shift approach for 
closed systems are two established definitions of conduc- 
tance, which compare rather favorably with each other 
for electrons in disordered systems (see Ref. 

The search for universal features in the statistics of sta- 
tionary electronic levels and their motion as a function 
of some parameter is a large area of study with obvi- 
ous consequences to quantum transport. For example, 
avoided level-crossings which characterize the spectrum 
of a chaotic system under a perturbation, such as the ap- 
plication of external fields, evolve with universal random 
matrix theory (RMT) laws j^. The very small level- 
spacings for avoided level-crossings correspond to large 
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absolute level-curvatures so that the tail of the level- 
curvature distribution P{K) is intimately related to the 
small part of the level-spacing distribution. Therefore, 
in the absence of symmetry breaking mechanisms the 
linear part of the Wigner distribution for small level- 
spacings implies large level-curvatures having a decreas- 
ing asymptotic behavior P{K) ~ 1/\K\^^^, with (3 — 1 
0^0, 1^ ■ This has been confirmed experimentally for 
level-curvature distributions of acoustic resonance spec- 
tra from quartz blocks . 

We focuse on the critical distribution of level- 
curvatures by obtaining the response of a quasiperiodic 
complex system's energy spectrum to different BC 10], 
within Thouless's approach to quantum transport. For 
this purpose a phase factor exp(i0) is imposed to the hop- 
ping probability connecting the first and the last sites of 
a ring and the resulting parametric dependence of the 
energy levels to changes of (j) is obtained. This allows to 
explore the nature of the corresponding eigenstates from 
extended to localized, since extended states simply feel 
any changes in BC having large curvatures while local- 
ized states are insensitive to BC having curvatures which 
approach zero. Nevertheless, the studied model is inter- 
esting only at the metal-insulator transition point since 
it satisfies duality for the extended fquasiballistic) and 
the localized (nonrandom) regimes One may won- 

der what the level-curvature distribution might be at the 
metal-insulator transition for such a quasiperiodic system 
which has Cantor set-like fractal electronic structure and 
a semi-Poisson hybrid distribution of level-spacings ^3 ■ 

Deviations around the maximum peak between level- 
curvatures obtained in experiment 9] and RMT were at- 
tributed to the presence of hidden symmetries 0, 0| . 
In our study the critical distribution of level-curvatures is 
obtained in a convenient one-dimensional setting which 
allows to discuss both the three-dimensional disordered 
system and questions of universality in general. Our find- 
ing is a scale-invariant P{K) which resembles the distri- 
bution obtained for the normalized level-curvatures in 
3-D critical disordered systems 0, having both a 
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diffusive tail and overall localized behavior For a 
finite system, the transition from extended to localized 
states involves broadening of the distribution P{K) and 
lowering its maximum peak as the curvatures move to 
lower values by increasing disorder. Our study, on one 
one hand, is in agreement with the semi-Poisson level- 
spacing distribution for the same system |lOl |. On the 
other hand, the coexistence of the localized almost log- 
normal broad form with the diffusive tail might be the 
reason for the lowering of the maximum peak observed 
in the experiment j^. 

We have studied the distribution function of level- 
curvatures, defined as 
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for the energy levels £q of a quasiperiodic ring in the pres- 
ence of phase exp(j0), obtained from the critical tight- 
binding model Hamiltonian |lO| 

n (nm) 

- X! (exp(i(?!))cJjC™ + exp(-i(/))c]'„c„) . (2) 

(nm) 

The sum is taken over all sites n, c„ (c^) is the anni- 
hilation (creation) operator on site n and the potential 
Vn = 2cos(2'Kan) is chosen at criticality with a the golden 
mean irrational. The eigenvalues Sa have correspond- 
ing eigenvectors \a) = J2x=i l V'a(2;)|a;) with amplitudes 
'4'a{x) = {x\a) for a finite chain of size L = Fi with 
(J = Fi-i/Fi, the ratio of two successive Fibonacci num- 
bers -Fi-i, Fi. The imposed general boundary condition 
to the wave function ^p{x + L) — e'^'^'ip{x) is equivalent to 
piercing the ring by the Aharonov-Bohm magnetic flux 
<i> via (f) — 27r<i>/<i>o, with $0 = ch/e the flux quantum. 

The curvatures Ka are properties of the energy levels 
in the limit of zero flux so that perturbation theory can 
be applied :4]. Using the small-0 expansion exp{±i(l)) = 
1 ± 10 — 0^/2 -I- ... the relevant terms up to 2nd order 
are < a\V\a >= -0^ < a\l >< l\a > and < a\V\f3 >= 
i(l){< a\l X N\I3 > - < a\N >< >) so that the 
level-curvature for the eigenvalues is 
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a — 1,2,...,L, in terms of the eigenvalues Sa and eigen- 
vector amplitudes ^pa of the unperturbed ring (0 = 0). 
The level-curvatures Ka give quantum transport prop- 
erties and they are more difhcult to obtain than level- 
spacings since one also needs the eigenvectors of the 
Hamiltonian. The formula of Eq. (3) is exact as long as 
2/3 7^ £a, since higher orders vanish at (f> = 0. Moreover, 



it is also both conceptually and numerically meaningful 
since the wave function amplitudes at the two ends of the 
chain V'all), 4'a{L) are expected to play a role for trans- 
port, while linear algebra relates the first term in Eq. (3) 
to the error in eigenvalue determination from truncating 
an infinite tridiagonal matrix . 

To use Eq. (3) we only need the eigensolutions of the 
periodic BC problem, that is the solutions of the Hamil- 
tonian of Eq. (2) with = 0. The problem simplifies 
further by the presence of the parity symmetric potential 
cos{x) = cos{—x) which allows to find symmetric and an- 
tisymmetric solutions separately, dealing with two tridi- 
agonal matrices instead of the full matrix in the presence 
of BC Tt^. This dramatic reduction in storage require- 
ments is achieved if the one-electron basis is rearranged 
to run within — s,— s-fl,...,s — l,s also by distinguishing 
between even and odd size L. For example, for odd L — 
2s -|- 1 the problem is reduced to a tridiagonal matrix of 
size s-\- \ for the symmetric states which has diagonal ma- 
trix elements the potential values Vqj ^s-i, Vs + 1) 
unity elements lying next to the diagonal, except in the 
first (second) row second(first) column where is \/2, and 
each eigenvector amplitude (x(0), x(l), x('S)) is related 
to the original via the symmetric rule V'(O) — x(0)i 
■0(— s) = V-'(s) = x{s)/ y/i^)- The second tridiagonal ma- 
trix for the antisymmetric states is of size s with matrix 
elements Vi, V^-i, — 1 on the diagonal and unity 
next to diagonal with corresponding eigenvector ampli- 
tudes (x(l), x{s)) related to the original via antisym- 
metry V(0) = 0, V'(-s) = -V'(s) = -x(s)/v/(2). By 
replacing with x Eq. (3) splits into a symmetric part 
labelled by a = 1, 2, s + \ 

i^. = (x.(.))^ + 2 5' Mf)x^, (4) 
and an antisymmetric part for a = .s + 2, s + 3, 2s + 1 
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In Eqs. (4), (5) the s-th last element of each eigenvector 
can be computed iteratively with no need to increase the 
storage requirements beyond that of a tridiagonal ma- 
trix. Moreover, the corresponding sums over (3 run over 
opposite kind of symmetric and antisymmetric states, re- 
spectively, being precisely zero for species of the same 
symmetry. A similar formula for even L is easily ob- 
tained. 

A ballistic ring has K^ = (Stt^/L^) cos(27rQ;/L) which 
give square-root singularities for P{K). These disap- 
pear for the critical ring studied where the curvature dis- 
tribution P{K) is much more complex so it turns out 
more convenient to consider the logarithmic distribution 
P{\a\K\) = \K\P{K), instead. The computed critical 
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FIG. 1: The computed critical distribution for the logarithm 
of level-curvatures P(ln li^l) and various sizes L is shown to 
be almost scale-invariant. The histograms for each size are 
obtained from the diagonalization of Eq. (2) employing the 
formula for the level-curvatures of Eqs. (4), (5). The continu- 
ous line is a Gaussian fit of the data. 

distribution presented in Fig. 1 for various sizes L is 
shown to be scale- invariant, overall approximated by a 
logarithmic normal form. In Fig. 2 we present a log-log 
plot of P(ln|iir|) which demonstrates the diffusive uni- 
versal tail P{K) ~ \K\-^ for large-|i4:|. For stcl&\\-\K\ 
we approximately find _P(ln|i^|) ^ \K\'^ . Unfortunately, 
rapid loss of accuracy for the too small curvatures does 
not permit to extract reliable exponents from fits of P{K) 
to the form suggested in Jjy . However, the overall behav- 
ior of the critical distribution shown in Fig. 3 is roughly 
similar to the critical distribution obtained for the ZD 
Anderson model 0, . 

The level-curvatures in disordered or chaotic systems 
exploit a sort of "dynamics" of the quantum stationary 
spectra as a function of flux which might be thought 
of "time". The conductance is then obtained from the 
level-curvatures in terms of eigenvalues and eigenvectors 
of the tight binding system with periodic BC. In our case 
the computational effort is minimized because the par- 
ity symmetry of the potential reduces the problem to two 
simple tridiagonal matrices for the determination of sym- 
metric and antisymmetric states. The needed last ele- 
ment of each eigenvector can be also computed efficiently 
with no increase of storage so that the size of the matri- 
ces can easily exceed 10^. We remind the reader that for 
diffusive disordered systems the computed averaged LB 
conductance {g) was shown to be related to the mean 
absolute curvature via {g) = TrdX |)/A, where A is the lo- 
cal mean-level spacing. For localized disordered systems 
the curvatures diminish and both distributions approach 
a log-normal form with (hig) = 7r(ln(|if |)). Our study 
of the unornialized level-curvatures for the quasiperiodic 
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FIG. 2: Log-log plot of the critical level-curvature distribution 
for various sizes L as in Fig. 1. Apart from the less accurate 
points in the left corner we see a slower behavior for small- 1 A' 
by plotting the equivalent of P{K) ~ \K\, while for large-l-Zf 
the diffusive asymptotic tail becomes P{K) ~ 1""^. 
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FIG. 3; Another demonstration of the critical level-curvature 
distribution for two sizes L — 17711, 28657 is shown which 
resembles the distribution for a 3-D disordered system [l^ll5|. 
The distributions for the quasiperiodic and the disordered 
system are rather similar although in the disordered case the 
curvatures are normalized divided by the local mean level- 
spacing A. 

system enabled us to obtain the critical P{K) for very 
large system sizes. The curvatures become very small for 
such large systems which set limits to the accuracy of the 
distribution. Our results are similar to that obtained in 
iD critical disordered system. Moreover, such computa- 
tions might suggest that the lowering of the maximum 
peak observed in experiment j^] could be thought of as 
due to an approach towards the critical region. Although 
the critical distribution shown in Fig. 3 still displays a 
diffusive tail of the log-normal form it is much broader 
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having less height than the pure diffusive RMT one. 

The obtained results for non-interacting fermions in 
a quasiperiodic ring complement and confirm previous 
studies at the metal-insulator transition of disordered 
systems. The distribution P{K) for the ensemble of crit- 
ical states at the transition (there is no ensemble over 
disorder) is scale-invariant like the corresponding level- 
spacing distribution which is known to depend sensitively 
on BC JilJ. After average over BC the critical distribu- 
tion of the level-spacings can be described by the semi- 
Poisson curve which combines both extended and local- 
ized behavior. Similarly, the obtained P{K) exhibits a 
hybrid character like what is obtained for critical 3D dis- 
ordered systems. In closing, it is remarkable that such a 
simple one-dimensional non-random model can capture 
most features displayed at the Anderson metal-insulator 
transition of realistic disordered systems. 
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